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Abstract 

We present a conjecture for expressing the coefficients in the Cayley- 
Hamilton theorem for supermatrices in terms of supertraces. The con- 
jecture is tested for several supermatrix dimensions and unique results 
are obtained. Generating functions for determining these coefficients 
are also given. 
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In the classical theory of matrices the Cayley-Hamilton theorem asserts that 
a given matrix M satisfies its characteristic polynomial, the determinant of 
the matrix {xI — M). For an n-dimensional matrix this is a monic polynomial 
of degree n in x, whose coefficients can be expressed as polynomials in the 
traces of the first n powers of the given matrix. The generalization of the 
Cayley-Hamilton theorem to supermatrices is not straightforward because, 
when M is a supermatrix, the superdeterminant of {xl — M) is a ratio of 
polynomials in the variable x. 

Let the (even) supermatrix M have the block form 

where A, B, C, D are, respectively, pxp, pxq,qxp, qxq dimensional 
matrices with p + q = n. The submatrices A, D have even and B, C have 
odd parity. We will denote by a{x) {d{x)) the characteristic polynomial of the 
submatrix A {D). Then the characteristic function (the superdeterminant of 
{xI—M)) can be written as the ratio of two polynomials in x in two equivalent 
ways P, [2] : 

h{x) = sdet{xl - M) 

det[d{x) {xl -A)-B adj {xl - D)C] 



det[a{x) {xl - D) - C adj {xl - A)By 

where adj(M) denotes the adjoint (=determinant times inverse) of M. 

In [2], the polynomial V{x) = a{x)'^~^^d{xy~^^ has been proposed as the 
characteristic polynomial, which is of degree 2pq + p + q in the variable x. 
When a{x) and d{x) have no common factors, it is shown in [1] that h{x) can 
be written as {a{x) + r)/ {d{x) + 1) where r, t are polynomials with zero body 
and degree less than a{x), d{x), respectively. In this case the characteristic 
polynomial becomes P{x) = {a{x) + r){d{x) + t), which is of degree p + q, 
i.e., the dimension of the supermatrix, just as for ordinary matrices. 

The definition h{x) = sdet{xl — M) is invariant under general similar- 
ity transformations M Um M U^l . When the similarity transformation 



d{x] 



a[x 
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supermatrix Um is block-diagonal, the polynomials appearing in the numera- 
tor and denominator of h{x) are separately invariant o This implies that the 
coefficients of x in these polynomials are invariants of the supermatrix un- 
der similarity transformations with a block-diagonal supermatrix Um- One 
expects that, as for ordinary matrices, these coefficients will be expressible 
in terms of the supertraces of powers of M: Sj = str{M^). If M satisfies 
identically a characteristic polynomial of degree Umax, then the supertrace of 
this identity implies that the supertraces of the ffist Umax powers of M must 
appear in the coefficients. Higher powers of M and corresponding supertraces 
will then be expressible in terms of these using the polynomial identity. But 
the supertraces are polynomials in the eigenvalues of M whose number equals 
the dimension of M , so rimax must be equal to p + q. 

Observe that the coefficients of x in the polynomials of the numerator 
and denominator of h{x) are themselves polynomials in the matrix elements 
(only determinants and cofactors appear - no inverses); and recall that the 
characteristic polynomial V{x) = a{x)'^~^^d{xy~^^ is of degree 2pq + p + q, 
while the reduced polynomial P{x) = (a(x) + r){d{x) + t) of degree p + g is 
obtained when a{x) and d{x) have no common factors. These considerations 
lead to the following conjecture: 

The Cayley-Hamilton Theorem for supermatrices: A (p, q) supermatrix 
M, whose diagonal submatrices (A, D) have eigenvalues (a^, di) with distinct 
bodies 7^ 5^), satisfies the equation 

a2p,M" + a2p,+iM"-i + .... + a^p^+J = 0, (4) 

where n = p + q and the are polynomials in the supertraces Si, Sn of M, 
which are homogeneous of total degree k in the elements of the matrix M, 
so that every term in the sum is of degree 2pq + p + q in the matrix elements. 

Using Mathematica, we have tested this conjecture for all values of (p, q) 
satisfying p + g < 6, as well as (5, 1) and (1, 5). In all cases, the supertrace 
polynomials ak{Si) are determined uniquely (up to an overall numerical fac- 
tor) by the assumptions stated. In fact, these polynomials can be determined 
using a diagonal {p, q) supermatrix; and subsequently it can be verified that 
any such (p, q) supermatrix satisfies (jlj) with the same expressions for as 

^ These statements follow from the definition of these polynomials in terms of deter- 
minants and the properties of determinants: det{xl ~ UAU^^) = det{U{xI — A)U^^) = 
det{xl — A), which are true also for superdeterminants. 
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functions of the S^. In all cases examined the coefficient 02^5 of the highest 
power of M is the square of a polynomial of order pq in the matrix elements, 
while the coefficient a2pq+i is the product of this polynomial with a second 
one (of order pq + 1). We give below the complete expressions of the super- 
Cayley-Hamilton identities for the lowest supermatrix dimensions {p+q < 5) 
choosing the coefficient of Si^^"^ M^+^ to equal unity: 

(1.1) : Si'^M'' -SiS2M + ^[-Si^ + S2'')l ^0, (5) 
(2, 1) : {S^^ - ^2)' M3 - 1 - S2) {S,' + 3 5i ^2 - 4 ^3) 

- ^ {Si^ - 9 Si^S2 + 16 Si^Ss - 9 Si'' 82^ + 9 ^2' - 8 ^3') M 

+ ^{S,^-AS,S3 + 3S2'){Si^~3S,S2 + 2Ss) 1 = 0, (6) 

(1.2) : [Si^ + ^2)' + ^ (^i^ + ^2) (^1' -SS1S2-4. S3) 

- 4 (-^i' + 9 S^^S^ + 16 -51=^53 - 9 Si''S2^ -9S2^-8 S3') M 

lo ^ ' 

- ^{Si'-ASiS3 + SS2''){Si' + 3SiS2 + 2S3) 1 = 0, (7) 
(3,1) : (5i-^-3 5i^2 + 2^3)'m^ 

- Usi''-^SiS2+ 2S3)(Si''-ASiS3-3S2' + 6S4)M^ 

+ ^ [-^i' + 4 5i*^^2 - 32 51^53 - 6 5/ (^2' - 6 54) 

+ 64 S^''S2 -53-4 S^' (9 S2^ + 18-52-54 + 8 ,53') 

+ 96,51-52^,53 + 9,52^- 32,52,53^ + 36,54^ -36 ,52^,54] 

+ ^ [^i^ - 12 ,5i^,52 + 24 ,5i«,53 + 18 ,51^ (,52^ - 2 S^)) 

+ 24^/52^3-12 51^(3 52^-6 52 54 + 8^3^) (8) 

- 72 51^53 (52^ - 2 54) - 8 53 (9 52^ - 18 52 54 + 8 53') 
+ 3 5i (27 52^ - 36 52^54 + 32 52 53^ - 36 54^) ] M 

- TTr (-^i^ - 6 '5i'52 + 3 52' + 8 5i 53 - 6 54) 

X [5i^ - 3 5i^52 + 9 51^52' + 9 52^ - 24 5i 52 53 
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- 8 ^i^^3 + 16 ^3 + 18 5'4(^/ - ^2)] / = 0, 
(2,2) : [Si^~AS,Ss + 3S2^y 

- 2(Si^-ASiS3 + 3S2'') (Si^S2-3SiS^ + 2S2S3]M' 



+ [-^S^'' + ls^'S2' + S^' (4 S3' - 9 S2 S,) - fs, S3' 

+ ^Si' (^2' + 4 54') + 3 ^2' (4 S3' - 3 52 S,)]M' 

+ [^S,^S2-S,'S, + S^(8S3'S,-^S2'-9S2S^'^ (9) 

+ 2 Si^S2 5*3 + 2 5i^53 ^4 + 5i^52 (s S2 S4 — 8 ^3^^ 

+ 6 Si'S2'S3 - l^i^Sa^ + 1^253 (9 ^2^4 - 8 S3')]M 



^[(5i^ + 9 51^52' - 8 51^53 + 16 53' - 18 52 ^4 
9 (5i^52 - 3 52^ + 8 5i 52 53 - 6 51^54)^] / = 0, 



We observe that the expressions for dimensions {q, p) can be obtained from 
those for (p, q) by changing the signs of all supertraces, Sj — )■ —Sj, as would 
be expected from the definition of supertrace. Moreover, as suggested by 
the fact that the coefficient of the identity always factors, we find that these 
matrix identities can be written as the product of two matrix polynomials 
of degrees p and q, respectively. This might be expected from the form of 
the reduced polynomial P{x) = (a(x) +r){d{x) + t), but it is not an entirely 
obvious conclusion. The (2, 1) identity, for example, is the product of the 
matrix polynomials P2,Pi: 

P2 = (5i2-52)m2-^(5i3-53)m+^(5i^-4 5i53 + 3 52') /, 
Pi = (5i2-52)m+^(5i3-35i52 + 253) /. 

The expressions for higher {p, q) values are exceedingly complex. How- 
ever, for supermatrices with particular symmetries it is possible that all co- 
efficients in (jl]) have a common factor - a polynomial in the supertraces - 
leading to an identity of total degree lower than 2pq+p + q. This happens for 
the OSp{p\q) supermatrices Z. They are graded antisymmetric and satisfy 
QZ + Z^Q = 0, where Q is the graded symmetric OSp metric and Z'^ is 



5 



supertranspose of Z. More explicitly, assuming covariant Z indices, 



v4 + A* = 0, D-D* = 0, 5 = C*, J + J* = 0. (11) 

Then the OS'pij^q) matrices M = QZ with mixed indices are graded an- 
tisymmetric and the supertraces of all odd powers of M vanish. Thus, for 
example, putting 5*1 = = 5*3 in the general (2, 2) case dHD we find an iden- 
tity of 8th degree in the matrix elements of 05*^(212), after dividing out a 
common factor 5*2^ : 

(2, 2)o5p : S2'M^~S2S,M^ + ^{-S2'' + AS^)l = 0. (12) 

Trying to determine the polynomials for the general (2, 4) case requires 
massive amounts of memory, as it requires manipulating 22nd degree poly- 
nomials in 6 variables; on the other hand, if we impose the extra conditions 
appropriate to the OSp{2\4:) case, we find there is no unique solution for the 
ttk in the super Cayley-Hamilton equation (jl]). However, the non-uniqueness 
is an overall factor of degree 8 in (jlj), resulting in a unique identity of 14th 
degree in the matrix elements of 05'p(2|4): 

2 - -6 



(2, 4)o5p : (^2' + 2 54) M' 



+ ^(^2' + 2^4) (^2' -6^2^4-16^6 
■[82^ + 18 32^3^-36 82^ S^- 
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72' 

+ 64 52^S'6- 72 128 S'e^jM^ (13) 

^ -[32' + 6 32^3^ + 12 32^3^'' + 72 ^2 ^4' 



144' 

8 32^36 - 96 32^ 3i + 96 3i^3e - 128 S.^r^I / = 0. 



When the A, D submatrices have common eigenvalues, the identities ob- 
tained above are still meaningful, although the coefficients of the different 
powers of M have zero bodies. Thus these identities cannot be used for 
expressing higher powers of M in terms of lower ones. It seems likely (see 
the examples in |2J) that for degenerate eigenvalues, identities with matrix 
powers higher than p + q may be needed, but never more than 2pq + p + q. 
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Having convinced ourselves of the validity of conjecture (jl]), we find gen- 
erating functions for calculating the coefficients. We present the results with- 
out proof since they were obtained by generalizing formulas obtained using 
diagonal supermatrices. We first write the characteristic equation (jl]) for 
{p, q) supermatrices as (assuming that the leading coefficient a2pq has non 
zero body) 

P+q 

bf'^^ MP+'^-^ = 0, bi"''^ = 1. (14) 

j=0 

The coefficients b^''^\S) are determined by a generating function F^^'''\S,t): 



1 r 



diF^P'^XS, t) 



j=0 J 

The generating function for the case q = (ordinary matrix) i^ [3] 

F(P'O) {S, t) = G{S, t) = e- S-i ^ . 
For g > we find that the generating function takes the form 



(15) 



(16) 



F(^'«)(5,t)= l-5:/i,tM G{S,t), 



k=l 



where the /i^'s are determined by the linear system of equations 









( bp+1 \ 




/ bp bp-i 


■ bp^qj^l \ 


B 






bp+2 


B^ 


bp+1 bp 


■ bp-q+2 




\f^qJ 




\ bp+q J 




\bp+q-l bpj^q^2 


■ bp ) 



(17) 



(18) 



In this equation bj is an abbreviation for b^j'~^'^'^\S). We have checked that 
the matrix B becomes singular only when an eigenvalue of submatrix A 
coincides with one of D, so generating functions can be obtained only for 
the non degenerate case assumed in writing (jlj). We also assume that q < p. 



The coefficients b^f''^^ for g > p are obtained from the corresponding bj^'^'{S) 
by reversing the signs of all Sj's as noted above. 



^ It is obtained most easily using the identity det{e^) 
and expanding the logarithm around a; = 0. 



„TrZ 



with Z = log(/ - xA) 
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The explicit form of the generating function for g = 1 is 



t) = (1 - ^1 tr G{S, t), = J (19) 

bp ' '{S) 

and for g = 2 is 

F'^P''\S, t) = (l-ii,t-ii2 G{S, t), (20) 

with 
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/,(p+2,0)x2 ,(p+2,0),(p+2,0) ' /,(p+2,0)x2 , (p+2,0), (p+2,0) " 

[Op ) — Op_i Op^i [Op ) — Op_i 

(21) 

It should be pointed out that in equations (IT71) to (PTi) all Sj's that appear 
in the final expressions for the generating function F^P''^\S,t) must be in- 
terpreted as supertraces str{M^) of the relevant [p, q) supermatrix M. For 
example, the Si in 63"^'°^ = 1(— S'l'^ + 3 6*2 — 26*3), which enters in both 
F^^'^\S,t) and F^'^''^\S,t), are different quantities in the two cases. 

These b^^''^\S) coefficients, obtained from the corresponding generating 
functions, are rational functions of the supertraces. The polynomial coef- 
ficients Ofc in (jll) are obtained by multiplying the corresponding b^''^^^s by 
(deti3)^, as the generating function is at most quadratic in the /ifc's. The 
first coefficient a2pq = (det in (jlj) is indeed homogeneous of order 2pq in 
the matrix elements, because the b^''''^^^s are homogeneous of order j. Simi- 
larly, the coefficient a2pq+i, depending on the term linear in t in the generating 
function, always has an overall factor of det B. 

The results reported in this note present strong evidence for the validity of 
conjecture (jl]). The generating functions obtained in flTTI) and (ITHi) are simple 
and can be used to derive general properties of the coefficients. It would be 
most desirable if this conjecture could be proved, and if the expressions for 
the coefficients given above could be derived rigorously. 



Addendum 

After this work was completed we became aware of similar results in the 
mathematical literature [H [5] . However, the final matrix identities obtained 
in these references, due to their derivation, contain the supertrace of one 
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higher power (S'p+^+i) than the dimension of the supermatrix. As a result 
these matrix identities, as given, cannot be factored into a product of two 
matrix polynomials of degrees p and q. Of course, when S'p+g+i is replaced 
by its expression in terms of 5*1, 5*^+^, the form of our identity ^ (up to 
a factor det B) is obtained. 
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